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We show that strong correlations play a vital role in protecting low energy excitations in disordered 
high temperature superconductors. The impurity-induced low-energy density of states (DOS) is 
greatly reduced in the strongly correlated superconductor compared to d-wave Bogoliubov-deGennes 
theory. The gapless nodal quasiparticles, and the resulting 'V in the low-energy DOS, are much more 
robust against disorder compared to the large-gap antinodal excitations. We discuss the relevance 
of our results to angle-resolved photoemission and scanning tunneling spectroscopy experiments. 



The effect of disorder on the strongly correlated d-wave 
superconducting state of high T c cuprates is a problem 
of great interest because disorder is intrinsic to doped 
Mott insulators. The resulting nanoscale inhomogeneity 
and its effects on the properties of cuprates has been 
investigated in a series of beautiful scanning tunneling 
spectroscopy (STS) experiments 0,0,0. 

One of the earliest puzzles regarding impurity ef- 
fects 0, in the cuprates was: why is d-wave super- 
conductivity (SC) so robust against disorder? It was re- 
alized over the years that there are at least two reasons 
for this: (i) the intrinsic disorder lies off the Cu02 planes 
in the spacer layers, and (ii) even in the case of inten- 
tional disorder, such as Zn on the planar Cu sites, the 
short coherence length in the cuprates leads to a very 
local and inhomogeneous response to impurities 0, in 
marked contrast to standard Abrikosov-Gorkov theory. 

In this paper we show that there is another important 
reason: disorder effects are greatly suppressed in the pres- 
ence of strong correlations, which are of course central to 
high Tc superconductivity. To make this point, we com- 
pare the results of two different T = calculations. (A) 
The first includes disorder-induced inhomogeneity in the 
superconductor using the Bogoliubov-deGennes (BdG) 
framework, but ignores correlation effects. In fact most 
of the theoretical work 0, H, S @ m tne field, barring 
a few exceptions uses just such an approach. (B) 
The second calculation includes the effects of both disor- 
der and correlations, combining the inhomogeneous BdG 
approach with the projection to states with no-double 
occupancy, as explained in greater detail below. By com- 
paring the results of (A) and (B) we find that: 

(1) Strong correlations lead to a renormalization of im- 
purity effects, resulting in a much shorter healing length 
over which the pairing amplitude is destroyed in the 
vicinity of an impurity. 

(2) The same impurity potential generates significantly 
fewer low energy excitations in the correlated system. 

(3) Gapless nodal quasiparticles are much less affected 



by disorder compared with the large-gap excitations near 
the Brillouin zone edge. 

(4) The robustness of the nodal quasiparticles can be seen 
by the very small impact of impurities on the "V" in the 
density of states N(w), i.e., on its \oj\ behavior at low 
energy, in the strongly correlated system. 

(5) The point nodes are protected in the correlated sys- 
tem in the presence of disorder. In contrast, in the ab- 
sence of correlations, the point nodes of the clean system 
expand into extended arc-like regions of gapless excita- 
tions induced by disorder. 

The robustness of the nodes and nodal quasiparti- 
cles as well as the nodal-antinodal dichotomy are of 
direct relevance to recent ang le-resolved photoemission 
spectroscopy (ARPES) fill Ha liH| and STS experiments 
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Model and Method: We use the t-t'-J model 
H = - } j t rT <{cl a c T : a + h.c.) 

r,r' ,ct 

+J (S r -S r ,-n r n r //4) + ^(F(r)-^)n r (1) 

<r,r'> r 

with a random potential V(r) to describe the effects of 
disorder in a strongly correlated superconductor. H acts 
on a Hilbert space with no doubly-occupied sites, re- 
sulting from the on-site Coulomb U, the largest energy 
scale in the problem. The kinetic energy describes elec- 
trons with spin a hopping from site r to r' on a 2D 
square lattice, with t the near-neighbor and t' the next- 
near-neighbor hopping amplitudes. The superexchange 
J = At 2 jU is responsible for d-wave pairing in the doped 
system and u fixes the average density n — 1 — x with hole 
doping x |l4j . The impurity potential V(r) is Vq > at 
randomly located sites with a concentration ni mp and is 
zero otherwise. We will focus on weak potentials Vq = t 
as a simple model of intrinsic disorder in cuprates. 

The most commonly used approximation is to sim- 
ply ignore the no-double-occupancy constraint and solve 
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the BdG mean field equations for an inhomogeneous SC 
0, ■ The interaction in eq. is decomposed into local 
bond pairing amplitudes A rr / = J(c r j<Vj — c r ^c r ^)/2, 
and Fock shifts W rr ' = J{cJ, CT c rCT ) /2 with the density 
n ( r ) = (So- c ro- c ro-)- The effective BdG Hamiltonian is 
then diagonalized and the fields defined above calculated 
self-consistently. The results obtained by ignoring corre- 
lations will be called 'BdG results' from now on. 

Our main goal is to include strong correlations over and 
above the inhomogeneous BdG theory. The Hilbert space 
consists of states $} = P|$o) where the projection oper- 
ator P = n r (l — n rT n rj.) suppresses double occupancy in 
any state |4>o)- We use here the Gutzwiller renormalized 
mean field theory fl5t Hi| to handle projection. In the 
translationally invariant, clean system such an approxi- 
mate treatment is in very good semi-quantitative agree- 
ment [To] ] with numerical studies of variational wavefunc- 
tions with P implemented exactly |ld|. 

We next generalize the Gutzwiller approximation to 
inhomogeneous states. The main idea is to write the ex- 
pectation value of any operator Q in a state P|$o) as the 
product of a Gutzwiller factor gQ times the expectation 
value in|J>o} so that (Q) ~ <7q(Q)o- The standard pro- 
cedure 16] for calculating gq can be generalized to keep 
track of the local density x(r) = 1 — n(r). We thus obtain 
the kinetic energy (cl a c r r a ) « <7t(r, r')(cJ cr c r ' (T )o and the 
spin correlation (S r -S r '} s» g s (r, r')(S r -S r /)o- We find the 
local Gutzwiller factors are given by g t (r, r') = g t (r)g t (r') 
with g t (r) = [2x(r)/(l + x(r))] 1 / 2 an d 5s ( r , r') = 4/[(l + 
x(r))(l+x(r'))}. 
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FIG. 1: Local density n(r) and pairing amplitude A(r) vari- 
ations, plotted on a normalized scale of zero to one, with 
Vo = l.Ot, riimp = 0.01 and average n = 0.8. Left panels: 
'BdG plus correlations' results (J = 0.33t). Right panels: 
'BdG' results (J = l.lt; see text). 

The BdG equations are then solved for 
eigenvalues E n eigenvectors (u„(r), v n (r)) together 
with local self-consistency equations for the density 
n(r) = 2 J2 n \v n (r)\ 2 , the pairing amplitude A r y — 
Mr, r') £>„(r>*(r) + u n (r)<(r')] and the Fock shift 
W r ,r> = J 2 (r,r')Y /n v n (r')v*(r), where Ji. 2 (r,r') = 
J(3 5s (r,r')±l)/4. 

Results which include correlation effects will be de- 



noted by 'BdG plus Correlations' below. To meaningfully 
compare results with and without correlations, we choose 
all parameters, except J, to be same: t = 1, t' = —0.25t, 
impurity potential Vq — t, and work at the same values 
of average density n and impurity concentration ni mp . J 
values are chosen such that we get the same A in the 
absence of disorder in the two calculations. 
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FIG. 2: Probability distributions (a) P(n) of density n(r), and 
(b) P(A) of pairing A(r) averaged over 20 disorder configu- 
rations for rii m p = 0.2 (other parameters same as in Fig. 1). 
no = 0.8 and Ao = 0.086 correspond to the case with no dis- 
order (ni m p = 0). (c) A(n) with and without correlations for 
nimp = 0, with J values chosen so that the corresponding A 
values match at n = 0.8. 



Local density and pairing amplitude: In Fig.^ wc 

compare the results obtained with (left panels) and with- 
out (right panels) correlation effects. In both cases dis- 
order makes the density inhomogeneous, with n(r) sup- 
pressed near the repulsive impurity. There is a marked 
difference in the spatial variation of the d-wave pairing 
amplitude A(r) = e rjr ' A r . r //4 where e r r > = 1 for 
r' = r±x and is —1 for r' = r±y. In the BdG calculations 
the impurity affects A(r) over a longer distance compared 
with the highly local response in the correlated system. 
One might think that this difference in healing len gths 
comes from the factor of two mass renormalization |lOl| 
due to correlations observed in the non-disordered case, 
however, we have checked that in fact the spatially 
varying Gutzwiller factors are responsible for the ob- 
served effect in Fig. 1. The suppressed electron density 
near an impurity modifies the local Gutzwiller factors 
and leads to an enhanced hopping on the bonds around 
it. The system is consequently able to repair the damage 
induced by the disorder and heal the pairing amplitude 
over a shorter length scale. 

We plot in Fig|2Ja) and (b) the normalized distribu- 
tions P(n) of the density n(r), and P(A) of the pairing 
A(r). P(n), which is a delta-function at no = 0.8 for the 
pure system, develops a bimodal structure in the disor- 
dered case. The low n peak comes from sites at or near 
the impurities, while the n > no peak leads to an average 
density Hq. We see that correlations produce a sharper 
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P(n) compared to the simple BdG calculation. In the 
Bdg calculation -P(A) is broadened over a range of A 
values strictly below the zero disorder A . In contrast, 
in the correlated system, P(A) has weight both above 
and below Aq. This difference in P(A) with and without 
correlations can be understood from A(n) variation in 
the uniform, clean system (Fig. Hfc)). The structure in 
P(A) in the correlated system can also be understood in 
detail, with various peaks associated with different kinds 
of sites: at, near and far from an impurity. 
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FIG. 3: DOS N(u>) for various n; mp [listed in panel (b)] aver- 
aged over 15 impurity configurations. Panel (a) shows results 
with and (b) without correlation effects. (Other parameters 
are the same as Fig. 1). The energy gap is 4Ao in the pure 
limit, (c) Low energy DOS integrated over a small window 



M < 0.02t as a function of in 



Results with correlations 



shown with open symbols (red) and plain BdG with closed 
symbols (blue). 



Tunneling density of states: In a disordered sys- 
tem the one-particle spectral function A(r, R; lo) — 
— ImG(ri, r 2 ; lo + i0 + )/ir depends on both the center- 
of-mass R = (ri + r 2 )/2 and the relative coordinate 
r = ri r 2 . Within the Gutzwiller approximation 
we obtain A(r, R;w) = g t (ri, r 2 )A°(r, R; u>) with A = 
£ n «n(ri)«*(r 2 )«(w-^) + «„(ri)«*(ra)«(w + ^). The 
Gutzwiller approximation describes the coherent part of 
the spectral function which dominates A for w's 
smaller than or comparable to the gap. We focus ex- 
clusively on low-energy excitations here [20|. 

STM experiments measure the local density of states 
(DOS) N(R,u) = A(0,R;w). We calculate the total 
DOS N(u>) obtained by averaging (N(R,u))r over the 
system and find that there are very significant differences 
between the results with (Fig.Qfa)) and without correla- 
tions (Fig. Ofb)). The most striking observation is that 
the behavior in the low energy DOS in panel (a) is 
hardly affected as n imp increases from 1% to 25%, in 
marked contrast with panel (b). Further, by integrating 
N(u>) over a small window \u>\ < 0.02£, we see in Fig.|3Tc) 
that very little low energy spectral weight is generated in 
the correlated system from the pair breaking effect of 
impurities compared with plain BdG calculation. This 
leads to one of our central results: low energy excitations 



are much more robust against disorder in the correlated 
system compared with the simple BdG calculation. 

Next we look at the DOS near the gap edge. In both 
Figs.GJa) and (b) the sharp log singularities ("coherence 
peaks") seen in the clean d-wave SC are suppressed with 
increasing n imp . The energy scale of the peaks in the 
DOS is shifted up slightly with increasing ni mp in the 
correlated case (see Fig. EHb)), while it seems more or 
less constant in the simple BdG results. 
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FIG. 4: Top: Low energy spectral function A(k,cu) for 
|w| < 0.024 for 200 x 200 system with n lmp = 0.2. The 
results with correlations (panel (a)) show only a slight ex- 
tension around the nodes, while the plain BdG results (panel 
(b)) show significant extension into "Fermi arcs", (c) A{k,uj) 
as a function of w for n; mp = 0.2 in the correlated system for 
two momenta: the node (in red) and the antinode (in blue). 



ARPES spectral function: Where do the the low- 
energy excitations seen in the DOS of Fig. [3] come from? 
To address this question we found it useful to com- 
pute the ARPES intensity, which in a disordered sys- 
tem is given by the spatially averaged spectral function: 
A(k, u) = E r expHk • r) (A(r, R; w)) R . In Fig. H we 
plot the near-zero energy intensity A(]s., ui) and see that 
the four point nodes of a clean d-wave SC are extended 
into "arcs" due to disorder ji^l- The gapless excitations 
at the k-points on these "arcs" contribute to the non- 
zero DOS. A similar effect has been previously found in 
a T-matrix calculation |22j which ignores correlations. 
The surprising new feature of our results is that the ex- 
tension of the "arcs" is much less when correlations are 
taken into account; compare Figs.^fa) and (b). This is 
yet another manifestation of the fact that disorder effects 
are suppressed in the presence of strong interactions. 

The calculated ARPES intensities also permits us to 
compare the behavior of the gapless nodal excitations 
with the maximum gap antinodal states. The coherent 
piece of the T = A(k, us) in a clean SC consists of 
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delta-functions at u> = i-Ek even in a strongly correlated 
system [l^. In the disordered SC these will be broad- 
ened by impurity scattering. We show in Fig. 0Jc) the 
^-dependence of A(k,uj) at two k-points: the node (N), 
along the zone diagonal, which is gapless and the antin- 
ode (AN), at the zone edge, which has the maximum of 
the d-wave gap. 

We see from in Fig. 0Jc) that the nodal spectrum is 
much sharper compared to the antinodal one. A simple 
golden rule estimate of the impurity scattering rate sug- 
gests that the linewidths should scale like the DOS at the 
excitation energy. Since we have already seen that the 
impurity-induced DOS in the correlated system is greatly 
suppressed relative the simple BdG result at low energies, 
it is self-consistent that the gapless nodal states are much 
less affected by disorder than the gapped antinodal ones. 

Our results on the nodal-antinodal dichotomy, ob- 
tained by including the effects of strong correlations 
and inhomogeneity, capture the essential features of 
both ARPES [H El E3 and STM @ experiments on 
cuprates. In STM experiments on BISCCO 0], a clear 
distinction is seen in the behavior of the local DOS for 
lo and for lu near the gap. The very low energy spec- 
tra are homogeneous and do not vary from one region 
of the sample to another. But the spectra near the gap 
are very heterogeneous, and thus sensitive to disorder. 
ARPES experiments [13, UM dearly see the greater sen- 
sitivity of the antinodal states to underdoping compared 
with the nodal states and this has been attributed to the 
appearance of competing order at low doping. Our re- 
sults indicate that disorder also has the same effect with 
antinodal excitations strongly affected and the nodal ones 
relatively robust. 

Conclusions: The calculations presented here focus on 
the interplay of strong correlations and disorder effects 
in superconductors, both of which are crucial elements of 
modeling high T c cuprates. We have shown that corre- 
lations suppress the effects of weak impurity scattering 
leading to a shorter healing length around an impurity, a 
reduction of the impurity-induced low-energy DOS and, 
most significantly, robust gapless quasiparticles. We have 
found that the nodes themselves are much less sensitive to 
disorder in a correlated system; in the absence of correla- 
tions the nodes expand into "arcs" of gapless excitations. 
The robust nodal quasiparticles lead to a 'V in the den- 
sity of states which is highly insensitive to disorder. We 
also found that the large gap antinodal excitations are 
much more affected by disorder than the gapless nodal 



ones. It would be very interesting to extend our results 
to other types of disorder - such as unitary scatterers 
and to hopping and superexchnage disorder Q - in the 
presence of strong correlations. 

We thank R. Sensarma for discussions. 
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